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Abstract 



We compute the gluon polarization AG and the gluon total angular momentum J g 
in the Isgur-Karl quark model. Taking into account self-interaction contributions, we 
obtain for both quantities a positive value of about 0.24 at the model scale /Xq — 0.25 
GeV 2 . This estimate is reasonably close to the expectation from polarized deep inelastic 
scattering analyses. 
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The EMC spin effect, i.e. the discovery that a surprisingly small fraction of the 
proton spin is carried by quarks, has been a major issue of hadronic physics for quite 
a long time by now Jl|. In the naive quark model one expects that quarks carry all 
proton's spin, that is AS = / dx (Aw + Ad + As) ~ 1. Taking into account relativistic 
effects this value is typically lowered to about 0.7 — 0.8, still much larger than the 
experimental finding for the singlet axial charge JJ. This discrepancy is what has 
been termed "the proton spin crisis". The solution of the puzzle is known 0. Due 
to the axial anomaly, the singlet axial charge, which is the experimentally measured 
quantity, is not simply AE but the combination ao = AS — (Nf a s /2ir) AG, where AG 
is the gluon helicity. A small value of ao (~ 0.3) is thus compatible with a relatively 
large value of AE if one allows for a sizeable AG. It has been shown J|] that the 
experimental data suggest AG = 1.6 ± 0.9 at Q 2 = 1 GeV 2 . The question arises 
whether a realistic model of the proton - or any non perturbative calculation - can 
predict a gluon polarization of such an order of magnitude. 

The situation is rather controversial. In || Jaffe gave a negative answer to the 
question: he found AG ~ —0.4 in the MIT bag model at a scale /x 2 , ~ 0.25 GeV 2 
and an even more negative value (~ —0.7) in the non-relativistic quark model. On the 
contrary, a recent QCD sum rule calculation || gives AG ~ 2.1 ±1.0 at 1 GeV 2 . Hence 
two different non perturbative approaches seem to give contrasting results. Another 
open issue [0] is the relative weight of AG and J g , the gluon total angular momentum, 
or - otherwise stated - the sign and the size of the gluon orbital angular momentum 

w 

The purpose of this Letter is to clarify these problems. We shall show that: i) the 
negative value for AG obtained by Jaffe is a consequence of neglecting "self-angular 
momentum" effects; ii) the quark model is not incompatible with a positive and sizeable 
gluon polarization; Hi) a specific calculation leads to a value of AG which is not far 
from what is required by the polarized DIS data and from what is obtained in the QCD 
sum rule computations; iv) J g is close to AG at the model scale, and its value is in 
good agreement with the finding of Ref. |7]]. 

In practical calculations we shall use the Isgur-Karl (IK) model ||, which is suffi- 
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cient for our purposes. We do not aim in fact at a precision calculation, but only at 
a qualitative study and an order-of-magnitude estimate. Moreover, the possibility to 
carry out analytically the calculations allows for a better control of the result. Our 
main conclusions would be unchanged should we use more sophisticated models of the 
nucleon. 

The angular momentum sum rule for the nucleon reads P, |Hj 

\ = \mQ 2 ) + L q (Q 2 ) + J g (Q 2 ) (1) 

where AS is the quark spin, L q the quark orbital angular momentum and J g is the 
total angular momentum carried by the glue. All quantities on the r.h.s. of eq. ([!]) 
are gauge invariant and depend in general on a scale /i 2 |TUf . We shall be concerned 
with the gluonic term in eq. ([I]), which is given by 

J g (^) = (p J | / d 3 r [fx (5(f) x 5(f))] 3 \p |) . (2) 

J g cannot be decomposed in a gauge invariant manner into an orbital (L g ) and a spin 
(AG) part written in terms of local fields. Nevertheless in the A + = gauge it is 
possible to give a local expression for the gluon in terms of the color fields || 

AG(/i 2 ) = (p t 1 / d 3 r 2 Tr j [5(f) x A(f)f + A ± (f) -S±(r)}bT>- (3) 

The renormalization of the operators appearing in eqs. (|2|, |3|) induces a scale de- 
pendence of J g and AG. However, in a quark model the matrix elements in eqs. (0, |^) 
do not depend on any scale. The model describes the nucleon at a fixed scale, which 
is usually very low (< 0.5 GeV 2 ). 

The color fields generated by the quarks are (a = 1, . . . , 8) 
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5(f) = ]T5 a (f) = -V£$>1; 

a=l a=l 



5(f) = J2 B a {r) = V x A a (r) 

a=l a=l 



where the potentials 



are written in terms of the single-quark charge densities 

Pi(r) = gS 3 (f-fi) 
and current densities (c = convective, m = magnetic) 
Uf = Ji(r)+jr(r) 



(4) 
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The expressions above are written in the non relativistic form, in view of their use in 
the IK model. 

We now make a gauge transformation to new potentials A 11 , enforcing the condition 
A + = 

A 0a (r) = $°(r) , 

A a (r) = A\r)-V [ Z d($ a (x,y,()- 

Jo 



Eqs. (|, D become 

Up 2 ) 
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EE/ d " r f x (p r I (%(f) x §*(f)) \ P t) 

i,j=l a=l 



(5) 
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AG(/i 2 ) = E E J d 3 r(P T I [EKf x ^(f)j + [jSftf) x V#(r* 

ij=l o=l I 



+ 



(6) 



where f a {r) = J Z d( <5> a (x, y, (). 

Upon integration by parts, the second term in eq. (||) vanishes, whereas the fourth 
term equals the first one plus a surface integral, which also vanishes in the model at 
hand. 

—* \ a —* . — * — * 

Setting Ef = -£-E\ and analogous relations for A®,Bf, and denoting by (-) c , (-)osf 
the color and orbital-spin-flavor expectation values, respectively, we get 

Up 2 ) = ~ E E(A-A, a ) c / d 3 r [fx (p T \E t (r) x B t (f)\p T) sf] 3 , (7) 

4 i,j=la=l 
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AGV) = \i2 E(A^) C / d 3 r( P u\2 [Mr) x A^f + A^f) ■ B ±j (r)\\p t) osf . 

4 i,j=la=l J I J 

(8) 

An explicit model calculation of eqs. (|7|, [8]) requires some approximation. For in- 
stance, one can replace the color fields by their expectation values in the nucleon 
eigenstate. In this case, if one includes the self-interaction terms with i — j, J g and 
AG turn out to be exactly zero. The reason for this is simple. Since the ground 
state of the nucleon is symmetric with respect to the exchange of any pair of quarks, 
the matrix elements in the integrals of eqs. ([7], §) do not depend on the particle in- 
dices i and j. Thus the integrals factorize out and J g and AG vanish exactly because 
J2i,j J2 a (K^P = 0. In || the self-interaction terms are neglected since they could 
diverge and therefore require to be renormalized. The result is a negative AG, which 
is a direct consequence of J2 a A® A| = —8/3 for 2 7^ j. 

In our calculation of eqs. (|7], |8|) we will include also the self-interaction contributions. 
Let us consider for defmiteness AG (the calculation of J g proceeds in a similar way). 
We insert in the matrix element of eq. (|8|) a complete set of intermediate states \/3), 
which include the orbital excitations of the three-quark system. In our non-relativistic 
calculation only positive energy states will appear. The matrix elements 

(0|^(f)|/?>, (P\Mr)\0) (9) 

for states \(3) with orbital mixed symmetry now carry a non trivial dependence on the 
particle index. As a consequence, the integral in eq. (|8|) no longer factorizes out and 
AG is not forced to vanish, even including the self-field terms. The procedure outlined 
above corresponds to taking into account static contributions to the self-spin and is 
equivalent to introducing a time-independent Green function for the quarks, expanded 
on the modes of the model Hamiltonian. 

Eq. (H) then reads (we omit the second term which turns out to be zero) 

AGV) 4EEE(\ S - A?) / d 3 r [<p T |4W> x ((3\A 3 (r)\p ]} osi f . (10) 

/ ij=l a =l (3 J 

Since the color electric field does not contain any spin-flavor operator, the intermediate 
states \f3) must have the same spin-flavor structure as the proton. Moreover, the totally 
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symmetric or antisymmetric intermediate states give again -independent matrix 
elements and therefore a vanishing contribution to eq. fllQl). 

Before showing our results, we need to discuss the convergence of the series obtained 
in the mode expansion. The latter is expected to be at worst logarithmically divergent, 
as it happens in QED for an isolated electron ||. Thus, in principle renormalization is 
required. On the other hand our calculation is performed in a non-relativistic effective 
model, which is supposed to be valid up to excitation energies of the order of few GeV 
at most. For larger energies relativistic corrections should be included. The situation 
here is rather similar to the one presented in the classical paper of Bethe on the Lamb 



shift |p.l| . In that case a natural cut-off was introduced and assumed to be equal 
to the electron mass. Higher energy contributions were left to be computed using a 
relativistic approach and renormalizing. Since the divergence was a logarithmical one, 
the precise value of the cutoff was irrelevant and a very good estimate of the energy 
shift was obtained. Here we adopt the same procedure, taking into account excited 
states having energies up to ~ 1 GeV above the ground state. 

We calculated analytically the contributions to AG and J g using the IK wave- 
functions [0]. The parameters of the model are fixed so as to reproduce the hadronic 
spectroscopy ||. The results are shown in Tab. 1, where the contributions are labelled 
by the energy quantum number M and by the orbital angular momentum £ of the 
excited states. They are also classified according to the symmetry of the wavefunctions 
with respect to the exchange of quark pairs. Notice that only the states with mixed 
symmetry give a non vanishing contribution, as we anticipated above. 

The resulting estimates for AG and J g are given by the sum of the contributions 
in the last two columns of Tab. 1. Only the first three energy levels are presented 
in the table. For AG we pushed the calculation up to the states with M = 3 and we 
found 0.03 a s . The apparently stronger convergence of the J g expansion can perhaps be 
related to its gauge-invariant expression. On the contrary, AG is in general non-local 
and the local expression here adopted is obtained choosing a particular gauge. It would 
be interesting to ascertain if, using directly the non-local expression which introduces 
a natural cut-off, the expansion do converge. 
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The values of J g and AG have been expressed so far in terms of the strong coupling 
constant a s , which is fixed in the IK model by reproducing the A — N mass splitting. 
This gives a s = 0.9 ||, which corresponds to a model scale — 0.25 GeV 2 . The 
dynamical quantities computed in the model are implicitly taken at this scale. 

We can now present our final estimate 

Thus, at the model scale, very little room (if any) is left for the orbital angular mo- 
mentum of gluons: L g (nl) ~ 0. 

In order to see what happens at a larger scale, we perform a leading-order (LO) 
QCD evolution. This is of course a dangerous thing to do when the initial scale is very 
low, as in our case. However, the evolution of J g is known only at LO |T0| , [12], and we 
content ourselves with a semiquantitative result. J g and AG evolve differently [Hj. At 
LO the quantity a s AG is constant, whereas J g increases slowly, its asymptotic value 
being 8/(16 + 3N f ) = 0.32. We find 

AG(5 GeV 2 ) = 0.70, J 9 (5 GeV 2 ) = 0.31 . (11) 

Hence, as Q 2 gets larger, the gluon orbital angular momentum L g becomes negative 
and increases in absolute value. This scenario is similar to the one suggested by the 
QCD sum rule calculation of Ref. JF|. Of course, a consistent check of the angular 
momentum sum rule, eq. (|IJ), requires AE to be calculated in the presence of polarized 
glue. This calculation is in progress and will be reported in a future work. 

In conclusion, let us summarize our results. We showed that the quark model does 
not imply a negative gluon helicity, which would contradict the indications coming from 
polarized DIS data. In a non relativistic model, taking into account self-interaction 
contributions, we obtained a value of AG which is positive and close to the range of 
values prescribed by the phenomenological analyses of data. Finally, we estimated J g 
and found that it is almost equal to AG at the model scale. 

It is a pleasure to thank Z. Berezhiani, L. Caneschi, R. Cenni and U. Tambini for 
many useful discussions. 
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Table 1: Contributions to AG [eq. (|10D 1 and to J g from different intermediate states 
\j3). These are classified according to their energy and orbital angular momentum 
numbers, and their symmetry under quark pair exchange. 
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